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摘 要 : HovpE RE ROT 7r US AR BU CET) TEE BR A TUBE PE [8] , AS SCL METRIC 模型 为 例 ,基于 地 
表 温度 (了 ) 与 叶 面积 指数 (ZL47) 相关 趋势 图 ,按照 “7, 频率 直方 图 占 比 前 5% 的 最 高 温度 区 LAT 频率 直方 图 占 比 后 
5% 的 裸 土 区 为 “ 干 点 ' LT. 频率 直方 图 占 比 后 5% 的 最 低温 度 区 LAT 频率 直方 图 占 比 前 5% 的 湿地 或 农田 为 “ 湿 
点 “的 提取 方案 ,综合 确定 研究 区 的 干 湿 限 ,并 通过 下 垫 面 涡 度 相关 系统 实测 数据 验证 了 模型 佑 算 精 度 。 经 检验 ， 


校准 


期 和 验证 期 的 区 域 ET 估算 控制 精度 平均 为 15% ,与 以 往 研 究 成 果 相 比 ,结果 相对 合理 ,提取 方案 在 一 定 程度 


上 可 为 METRIC 模型 主观 确定 干 湿 限 提供 帮助 。 将 ET 估算 结果 按 土 地 利用 分 类 提取 ,研究 区 不 同 土地 利用 类 型 
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蒸 散 发 (evapotranspiration , ET ) 作为 陆 面 水 文 循 
环 和 能 量 循环 过 程 的 重要 环节 ,是 水 热平衡 联系 的 
纽带 ” ,准确 监测 和 估算 其 时 空 变化 规律 ,对 提高 
区 域 水 分 生产 利用 效率 、 缓 解 地 区 用 水 矛盾 等 具有 
重要 理论 价值 和 现实 意义 ””。 随 着 ET 研究 理论 
的 不 断 深入 ,涌现 出 许多 ,如 水 量 平衡 法 ARB DUI 
波 文 比 法 、 空 气动 力学 法 ` 涡 度 相关 法 闪烁 通 量 法 
等 监测 和 估算 方法 ,这 些 方法 在 微观 和 农田 尺 
度 可 详尽 描述 ET 过程 的 特征 规律 ,为 下 垫 面 水 汽 
交换 ,农田 灌溉 决策 水 资源 综合 管理 提供 了 重要 参 
HARI. ART, ET 研究 尺度 逐渐 由 微观 拓 
展 到 宏观 ,但 这 些 传统 的 监测 手段 和 方法 并 不 能 详 
尽 反 映 区 域 ET 的 真实 空间 特征 。 到 感 技术 的 出 
现 ,克服 了 传统 方法 中 定点 观测 难以 推广 到 大 尺度 
的 瓶颈 ,为 大 范围 ET 估算 提供 了 有 效 解决 方案 。 

近年 来 ,遥感 ET 模型 的 发 展 从 简单 到 复杂 、 从 
经 验 到 机 理 ,在 理论 创新 和 方法 应 用 上 已 取得 众多 
成 果 , 其 中 ,以 遥感 能 量 平衡 ( RSEB) 为 理论 基础 的 
METRIC 模型 是 当前 估算 区 域 ET 的 典型 代表 。 该 
模型 计算 的 关键 是 地 表 阻 抗 (7,) 和 地 气温 差 (d7) 
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两 个 参数 ,参数 r, 通过 莫 宁 - AAR AEE 
现 最 优 求 解 ; 另 一 参数 dT 在 求解 过 程 中 , METRIC 
模型 巧妙 利用 “ 干 点 ”和 “ 湿 点 ”两 个 干 湿 限 解 决 了 
经 验 参数 c 和 5 求解 问题 ,大 大 简化 了 显 热 通 量 计 
算 。 但 是 ,模型 要 求 研究 者 必须 对 研究 区 物理 背景 
进行 深入 了 解 ,确保 主观 选取 的 干 湿 点 为 研究 区 内 
显 热 通 量 和 潜 热 通 量 的 上 下 临界 点 ,才能 保证 显 热 
通 量 和 潜 热 通 量 的 计算 精度 。 人 研究 者 在 应 用 MET- 
RIC 这 类 模型 时 ,由 于 干 湿 限 的 确定 具有 很 大 的 主 
观 性 ,往往 因 找 不 到 理想 的 “ 干 点 ”和 “ 湿 点 ” 像 元 导 
致 模型 计算 精度 不 高 ,从 而 产生 遥感 ET 空间 歧义 
性 问题 和 25) , 即 不 同 研究 者 利用 相同 的 遥感 数据 和 模 
型 却 得 到 不 同 的 ET 估算 结果 。 其 原因 是 不 同 研 究 
者 认识 标准 不 一 致 ,主观 选取 的 干 湿 限 标准 不 一 , 计 
算 的 显 热 通 量 或 潜 热 通 量 差 别 较 大 ,从 而 影响 了 区 
Ja ET 模拟 精度 ,限制 了 模型 更 大 范围 的 应 用 。 基 
于 此 ,本 文 针 对 METRIC 模型 计算 区 域 ET 出 现 的 
空间 歧义 性 问题 ,通过 地 表 温 度 (7.,) 和 叶 面 积 指数 
(1L47) 梯 形 特征 空间 模型 理论 构建 干 湿 限 ,尝试 提 
出 一 种 “ 干 点 ”和 “ 湿 点 ” 像 元 的 简易 提取 方法 ,并 借 


究 院 科研 专项 


作者 简介 : 王 军 (1987 - ) , 男 , 工 程 师 ,博士 生 , 研 究 方向 为 牧区 水 资源 高 效 利用 与 流域 水 循环 模拟 . E-mail; slwj1988@ 163. com 


http ://azr. xjegi. com 


TOS 区 


助 下 垫 面 涡 度 相关 系统 实测 水 汽 通 量 数据 验证 模型 
估算 精度 ,进一步 探讨 ET 的 时 空 分 布 特征 ,以 期 为 
posi ET 精准 估算 提供 帮助 。 


1 研究 方法 


METRIC ( mapping evapotranspiration at high res- 
olution with internalized calibration ) 模型 是 以 SEBAL 
模型 为 理论 基础 17 ,考虑 下 热 面 高 程 .坡度 , 坡 向 
等 因素 的 干扰 ,通过 将 遥感 表面 信息 转换 为 地 表 反 
照 率 .地表 比 辐射 率 、 植 被 指数 和 T, 等 地 表 特 征 参 
数 ,结合 下 垫 面 气象 观测 数据 计算 区 域 ET 的 一 种 
Jk, He EA: 

LE-R,-G-H (1) 
式 中 :LE 为 潜 热 通 量 ;RR, 为 地 表 净 辐射 量 ;C 为 土壤 
热 通 量 ;五 为 显 热 通 量 。 式 中 各 能 量 分 项 单位 均 为 
W*m^, 
1.1 地 表 净 辐射 量 

R, 是 由 天 空 (包括 太阳 和 大 气 ) 向 下 投射 与 由 
地 表 ( 包 括 土壤 植物 水面) 向 上 投射 之 间 全 波段 
辐射 量 之 差 ,该 值 是 下 垫 面 能 量 水 分 传输 的 主要 动 
力 。 其 计算 方程 式 为 : 

R,=(-a) -K,, +(L, -£,,) -(-e) +L, (2) 


式 中 :a 为 地 表 反 照 率 ,无 量 纲 ;K;, 为 短波 辐射 ;L， 
为 下 行 的 长 波 辐射 ;L,, 为 上 行 的 长 波 辐射 ;se 为 地 
表 比 辐射 率 ,无 量 纲 。 式 中 能 量 分 项 单位 均 为 W， 
m”“。 各 参数 计算 依据 详 见 文 献 [14]。 
1.2 土壤 热 通 量 
G 是 指 由 于 传导 作用 而 存储 在 植被 和 土壤 中 的 

那 部 分 能 量 ,其 变化 主要 受 地 表 净 辐射 和 下 垫 面 植 
被 覆盖 影响 。 综 合 比较 ,本 文 参照 Bastiaanssen 和 
Wright ERRO ,采用 下 列 经 验 公式 进行 计算 : 

T. -273.15 
> Xs 


G (0. 003 8a +0. 007 4a’) - 


(1-0.978NDVI') - R, (3) 
IP:T, 为 地 表 温 度 (K) ;NDVI 为 归 一 化 植被 指数 ， 
其 计算 依据 卫星 热 红 外 波段 和 红外 波段 计算 ,无 量 
纲 ; 其 他 参数 意义 同上 。 
1.3 HS 
H 是 描述 陆 面 与 大 气 能 量 交换 的 一 个 过 程 , 通 
过 对 流 或 传导 的 方式 ,从 地 球 表面 转移 到 大 气 中 的 
那 部 分 能 量 。 对 于 五 的 计算 ,METRIC 模型 基于 “ 干 
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点 ”和 ” 湿 点 " 像 元 、 零 平 面 位 移 以 上 高 度 2 和 2 处 
温差 d7 (dT - T, - T) 5 T, Z RIZ TEORS AR UR, 
巧妙 地 避免 了 下 执 面 气象 数据 空间 插值 及 订正 7 
引起 的 误差 ,进而 通过 高 程 .坡度 , 坡 向 修正 ,实现 于 
感 影像 每 个 像 元 计算 结果 的 校准 。 其 计算 方程 式 
为 : 


dT 


H=p,,,¢, T (4) 


式 中 :pu 为 空气 密度 (kg * m 2) 50, 为 空气 定 奈 比 
热 容 (1 004 J+ kg! - KD ir, 为 两 个 近 表面 高 度 
之 间 的 空气 动力 学 阻抗 ,其 计算 是 假设 地 面 200 m 
以 上 的 摊 混 层 空 气 运动 不 受 下 热 面 影响 ,整个 研究 
区 域 为 同一 常数 ,利用 200 m 处 的 风速 推算 地 表 的 
摩擦 风速 求解 aT 通过 与 T, 之 间 线 性 关系 式 
计算 。 其 计算 方程 式 如 下 : 


dT x GT, adum 十 b (5 ) 
df. df 
a= : (6) 
T. adum? a T. dimi 
b = dT = aT, unt (7 ) 
enn = T. + 0. 006 5AZ (8) 


AP iab 为 经 验 参 数 ;Ti 为 修正 地 表 温 度 (K)， 
对 于 非 平 原 区 需要 利用 数值 高 程 差 AZ 修正 ,其 中 
AZ 为 每 个 像 元 的 高 程 与 参考 高 程 之 差 ;其 他 参数 意 
义 同 上 。 


2 研究 区 概况 与 数据 预 处 理 


2.1 研究 区 概况 

研究 区 选 在 内 蒙古 典型 草原 区 锡 林 河 流域 (地 
理 坐 标 为 115.53° ~117.25°E,43.41° ~44.63°N)， 
面积 约 为 11 172 km 。 地 势 呈 东南 高 西北 低 ,海拔 
在 902 ~1 621 m。 多 年 平均 降水 量 为 188 mm ,多 年 
平均 营 发 量 约 为 1 903 mm ,在 空间 分 布 上 由 东南 向 
西北 递增 , 属 典型 的 大 陆 性 温带 半 干 旱 气 候 。 区 内 
土地 利用 类 型 主要 是 以 草地 为 主 ( 图 1a 和 图 1b)。 
2.2 数据 源 及 预 处 理 

数据 源 处 理 包 括 遥 感 和 气象 数据 。 其 中 ,遥感 
数据 来 自 USGS 网 站 的 Landsat 5 TM 数据 (分 辩 率 
30 m) ,结合 云 履 盖 状 况 ,选用 2006 年 8 月 4 日 卫星 
过 境 时 刻 的 遥感 影像 作为 校准 期 数据 ,2007 年 7 H 
6 日 .2008 年 8 月 25 日 .2009 年 6 月 25 日 和 8 月 
12 日 作为 验证 期 数据 ,利用 DEM 对 原始 数据 进行 
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(a) 流域 所 在 位 置 
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(c) 涡 度 相关 系统 


图 1 


Fig.1 Geographical location of the study area 


表 1 Landsat 5 TM 数据 信息 
(PATH 124/ROW 29 和 30) 
Tab.1 [Information of landsat 5 TM Data 
(PATH 124/ROW 29 and 30) 


日 期 WREE TENE sAm 
/d 当地 时 刻 /% 
校准 期 ”2006 -08 -04 216 10:51:48 0 
验证 期 ^ 2007 -07 -06 187 10:52:15 0 
2008 -08 -25 238 10:44:09 0 
2009 -06 -25 176 10:47:00 6 
2009 -08 -12 224 10:48:09 0 


表 2 卫星 过 境 时 刻 各 站 点 气象 信息 
Tab.2 Meteorological information of the stations 
at satellite transit time 


经 度 ”纬度 海拔 风速 ”气温 水 汽 压 
£C) o) /m  /(m:s^) /CC /kPa 


站 点 


东 乌 珠 穆 沁 旗 116.97 45.52 840.50 2.8 25.5 0.91 
Bay ME 114.95 44.02 1 127.90 3.2. 23.9 0.89 
化 德 县 114.00 41.90 1 140.00 3.1 18.1 0.85 
Vi SERA ME 117.60 44.58 996.60 4.9 22.2 0.90 
锡林浩特 市 116.07 43.95 990.80 4.3 24.2 0.81 
林 西 县 118.07 43.60 800.30 0.9 23.0 0.92 


多 伦 县 116.47 42.18 1 245.40 1.3 19.3 0.88 


研究 区 位 置 图 


几何 校正 和 掩 膜 处 理 , 得 到 模型 所 需 的 遥感 数据 
( 表 1)。 

气象 数据 来 自 锡林浩特 市 等 周边 7 个 旗 县 的 气 
象 站 ,具体 包括 气温 、 风 速 和 水 汽 压 等 。 其 中 ,气温 
用 于 T, 和 地 表 净 辐射 计算 ;风速 用 于 空气 动力 学 阻 
抗 计算 ;水 汽 压 用 于 大 气 透 过 率 、 地 表 反 照 率 等 参数 
计算 。 以 2006 年 8 月 4 日 为 例 ,各 气象 站 点 在 卫星 
过 境 时 刻 (11:00 左右 ) 气 象 数据 整 理 见 表 2。 


3 PERSAS ET 估算 


3.1 干 湿 限 确定 与 分 析 

模型 中 的 干 湿 限 即 “ 干 点 ”和 “ 湿 点 ”两 种 特殊 
的 极限 像 元 。 其 中 ,“ 干 点 ” 指 没有 植被 覆盖 的 干燥 
的 闲置 农田 或 荒地 ,温度 很 高 , 蒸 散 量 几乎 为 0 的 像 
To, “Fa HR, -CE=0; 湿 点 "是 指 影 
像 中 水 分 供应 充足 植被 生长 茂盛 .温度 很 低 、 处 于 
潜在 燕 散 水 平 的 像 元 ,可 以 是 植物 生长 良好 的 完全 
窗 盖 的 区 域 或 开放 的 水 体 ,在 “ 湿 点 ”上 ~R, - 
G - LE ,特殊 情况 下 且 ~0。 通 过 对 “ 干 点 ”和 “ 湿 点 ” 
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点 位 信息 提取 , 便 可 以 求 得 dT 值 的 分 布 ,但 这 种 
提取 过 程 往往 伴随 空间 歧义 性 的 出 现 。 针 对 模型 计 
算出 现 的 空间 歧义 性 问题 , Tasumi 7? 借鉴 Moran 提 
出 的 dT7ZS47(S4 了 表示 土壤 调节 植被 指数 ) 梯形 
特征 空间 模型 理论 ,假定 dT GP ST T. 高度 相关 
的 前 提 下 ,参照 7. 与 LA 相关 趋势 图 选 定 干 湿 
限 (图 2)。 

图 2 中 ,A 点 附近 为 干燥 裸 土 区 ,集中 在 该 区 域 
的 像 元 表现 为 土地 裸露 缺 水 ;B 点 附近 为 湿润 裸 土 
区 ,集中 在 该 区 域 的 像 元 一 般 为 盐碱地 或 裸露 的 沙 


叶 面 积 指 数 (LAD 
2 地表 温 度 (7.,) 与 叶 面 积 指数 (147) 相关 趋势 图 
Fig.2 Correlation of surface temperature ( T,) and leaf 
area index (LATI) 
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地 ;C 点 附近 为 高 覆盖 湿润 区 ,集中 在 该 区 域 的 像 元 
一 般 为 水 分 湿润 的 湿地 或 农田 ;D 点 附近 为 高 覆盖 
BRAK X ,集中 在 该 区 域 的 像 元 一 般 为 植被 高 覆盖 的 
时 作 耕 地 或 草地 。 根 据 Tasumi 人 研究 结论 ,“ 干 点 ”在 
T. 频率 直方 图 占 比 前 5% 像 元 中 提取 ,“ 湿 点 "在 7 
频率 直方 图 占 比 后 3% ~ 5% 像 元 中 提取 。 

在 实际 操作 过 程 中 , 受 季 节 、 土 地 利用 类 型 及 云 
和 覆盖 等 因素 影响 , 仅 遂 过 7, 选取 干 湿 限 仍 存在 较 大 
的 偶然 性 和 不 确定 性 。 因 此 ,本 文 提取 方案 是 综合 
考虑 下 垫 面 T, 和 植被 覆盖 及 土地 利用 现状 ,首先 吻 
除 云 覆盖 的 干扰 ,在 Tasumi 提取 方案 基础 上 ,加 入 
L417 约束 ,将 “ 干 点 ”提取 范围 集中 到 A 点 ,“ 湿 点 ” 
提取 范围 集中 到 C 点 。 具 体 提 取 是 以 土地 利用 为 
底 图 ,考虑 有 效 像 元 提取 的 便利 性 ,LA 控制 线 设置 
过 小 ,符合 约束 条 件 的 像 元 难以 发 现 ,ZL47 控制 线 过 
大 ,符合 约束 条 件 的 像 元 过 多 ,无 法 提取 到 有 效 干 湿 
限 ,因此 LAT 干 湿 限 控制 线 参 照 工 中 的 5% 提 取 , 即 
选取 T, 频率 直方 图 占 比 前 5% 的 最 高 温度 区 LAT 
频率 直方 图 占 比 后 5% 的 裸 土 区 为 “ 干 点 ”, 选 取 T, 
频率 直方 图 占 比 后 5% 的 最 低温 度 区 LAT 频率 直方 
图 占 比 前 5% 的 湿地 或 水 田 为 “ 湿 点 ”。 
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图 3 研究 区 地 表 温 度 和 叶 面 积 指数 空间 分 布 


Fig.3 Spatial distribution of surface temperature (7,) and leaf area index ( LAI) in the study area 


图 3 是 区 域 “ 干 点 "和 “ 湿 点 " 像 元 提取 需要 计 
HAI T, ALALE. HEP, T, 依据 单 窗 算法 计算 ， 
经 式 (8 ) 高 程 调整 计算 ,2006 年 8 月 4 日 研究 区 
T_e [283.0 K,325.0 K ] ,均值 为 306.2 K (图 3a); 
LAI 依据 土壤 调节 植被 指数 (54V1) 计算" ,确定 研 
究 区 LAI e [0.0,6.0] ,均值 为 0.38 (图 3b)。 方案 
Adae A EE. T, 和 LAL 指数 图 Sms 


的 无 效 像 元 (图 3a P T, Æ 283 ~ 293 K 的 部 分 像 
元 ) ,确定 了 研究 区 的 “ 干 点 ” (T, =315.8 K,1AI = 
0.00) 和 “ 湿 点 ”(7T 2296.6 ,L141=1.73) 像 元 。 
3.2 结果 验证 

干 湿 限 确定 后 ,利用 莫 宁 - 奥 布 霍 夫 循环 迭代 ， 
实现 显 热 通 量 和 潜 热 通 量 计 算 ,进而 通过 蒸发 比 法 
PRE A ET 值 。 通 过 METRIC 模型 反 演 的 区 域 


王 “ 军 等 :基于 地 表 温 度 和 叶 男 


ET 是否 真实 地 反映 了 下 垫 面 蒸 散 耗 水 特征 ,需要 进 
行 结果 验证 及 对 比分 析 。 受 监测 技术 的 限制 ,目前 
仍 未 有 较 成 熟 的 技术 实现 ET 的 区 域 尺度 验证 。 对 
于 模型 的 估算 结果 ,本 文采 用 内 蒙古 锡林郭勒 草原 
生态 系统 国家 野外 科学 观测 研究 站 实测 数据 完成 。 
该 站 位 于 研究 区 东南 部 ,站 内 布设 的 涡 度 相关 系统 
实测 水 汽 通 量 数据 来 验证 METRIC 模型 的 估算 结 
果 。 经 现场 调查 , 涡 度 相关 系统 方圆 2 km 内 均 为 天 
然 草地 (图 16) ,植被 类 型 主要 以 大 针 茅 和 羊 草 为 
主 中 ,因此 可 以 根据 研究 区 土地 利用 现状 图 及 
NDVI 的 取 值 特征 ,选择 以 系统 所 在 位 置 (地 理 坐 标 
H 116. 67°E,43.55°N; Fik 1 250 m) 为 中 心 的 3 x 
3 大 小 窗口 范围 内 9 个 像 元 (0. 81 x 10* m^) ET fih 
算 结 果 平 均值 为 代表 ,进行 实测 的 蒸 散 量 与 MET- 
RIC 模型 估算 的 ET 之 间 的 对 比分 析 ( 图 4) 。 

图 4a 显示 了 研究 区 5 个 不 同时 相 的 METRIC 
模型 遥感 估算 ET 值 与 涡 度 相关 实测 ET 值 之 间 的 
数值 关系 ,2006 年 8 月 4 日 .2007 年 7 月 6 日 .2008 
年 8 月 25 日 2009 年 6 月 25 日 和 8 月 12 日 等 5 个 
不 同时 相 遥 感 估算 ET 值 与 涡 度 相关 实测 ET 值 的 
相对 误差 分 别 为 19% 27% 10% .4% 和 16% ,平均 
相对 误差 为 15% 。 结 合 Kalma 等 对 不 同 遥 感 估算 
ET 方法 的 精度 统计 发 现 ,与 基于 地 面 的 ET 测量 值 
相 比 ,大 多 数 相对 误差 范围 为 15.0% ~30. 0% 09 , 
从 误差 结果 分 析 , 这 种 干 湿 限 提取 方案 计算 的 区 域 
ET 估算 可 以 满足 区 域 ET 研究 的 精度 要 求 。 由 图 
4b 可 以 看 出 ,利用 METRIC 模型 估算 的 E7 与 涡 度 
相关 系统 实测 水 汽 通 量 数据 在 一 定 程度 上 变化 趋势 
整体 吻合 ,模拟 值 与 实测 值 之 间 具 有 较 高 的 拟 合 优 
度 (决定 系数 尼 为 0.98) ;根据 统计 显示 ,在 各 个 时 
相 上 METRIC 模型 估算 结果 均 存 在 不 同 程度 的 高 估 
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现象 。 受 仪器 设备 和 资料 数据 的 限制 ,如 果 需 要 对 
ET 估算 结果 进行 系统 的 、 多 区 域 的 对 比 验 证 , 则 需 
要 在 人 研究 区 布设 大 孔径 闪烁 仪 等 设备 ,来 进一步 印 
证 模型 估算 结果 的 可 信 度 。 
3.3 ET 空间 特征 分 析 

根据 METRIC 模型 的 基本 原理 及 算法 ,本 文 基 
于 自主 开发 的 区 域 ET 计算 系统 V1.0 (软件 著作 权 
登记 号 2017SR680304) ,计算 得 到 研究 区 不 同时 相 
的 ET7 值 (图 5)。 其 中 :图 5(a ~e) 依 此 为 锡 林 河 流 
域 2006 年 8 月 4 日 .2007 年 7 月 6 日 .2008 年 8 月 
25 H 2009 4Æ 6 H 25 H418 H 12 HH ET 结果; 图 
5f 为 研究 区 土地 利用 空间 分 布 。 

结合 研究 区 LA1( 图 3b) 和 土地 利用 空间 分 布 图 
(图 5f) 分 析 , 从 锡 林 河 流域 东南 部 的 林 草 交错 带 到 
西北 部 的 低 覆 盖 草 地 ,降水 量 逐 渐 减 少 . 干 旱 化 程度 
逐渐 加 剧 ,ZL417 也 存在 显著 差异 ,因此 在 气候 条 件 、 
土壤 供水 状态 、 土 地 利用 类 型 等 因素 的 综合 影响 下 ， 
锡 林 河 流域 ET 的 估算 结果 也 表现 出 明显 的 空间 分 
异 特性 。 将 图 5 中 ET 结果 与 土地 利用 类 型 相 结合 
进行 分 析 ,发 现 除 水 体 等 高 蒸 散 特性 地 物 之 外 ,5 个 
时 相 的 ET 估算 结果 均 表 现 出 了 从 东南 向 西北 递减 
的 变化 趋势 。 流 域 东南 部 为 大 兴安 岭南 夯 灌 木林 
区 ,降水 相对 丰富 , 林 草 植被 长 势 良好 ,具备 了 较 好 
的 有 7 条 件 。 另 外 ,近年 来 流域 内 少量 天 然 草地 开 
发 变 成 了 种 植 青贮 玉米 的 灌溉 人 工 草地 ,相对 充足 
的 供水 条 件 保 证 了 东南 部 林 草 交错 区 和 耕地 区 的 
ET 值 最 接近 于 水 体 的 蒸 散 强度 。 流 域 西北 部 属于 
典型 草原 区 ,该 区 降水 稀少 .气候 干旱 ,生态 环境 十 
分 脆弱 ,加 上 近年 来 过 度 放 牧 导 致 草原 类 型 逐步 退 
化 成 荒漠 草原 ,降水 稀少 和 植被 类 型 单一 共同 限制 
T 57, 因此 该 地 区 实际 ET 明显 低 于 其 他 区 域 。 
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Fig.4 Comparison between ET values estimated by METRIC model and measured by eddy covariance system 
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Fig.5 Estimated results of ET 


另外 ,将 5 个 时 相 的 ET 结果 平均 化 处 理 , 不 同 
土地 利用 类 型 表现 为 ET > ET yyy > ET yy > 


频率 直方 图 占 比 前 5% 的 湿地 或 水 田 为 " 湿 点 ”的 
提取 方案 ,综合 确定 了 人 研究 区 的 “ 干 点 ”和 “ 湿 点 ” 像 


Tj 发 和 用 t+ 地 > ET ge comin 于 Twig。 分 析 原 因 ， 
充足 的 水 分 为 水 体 ET 创造 了 条 件 ;林地 多 分 布 在 
河谷 低洼 地 带 , 下 垫 面相 比 其 他 类 型 水 分 较为 充足 ; 
耕地 多 是 在 天 然 草地 基础 上 开垦 形成 ,喷灌 和 滴灌 
等 灌水 技术 为 ET 提供 了 水 分 保证 。 另 外 ,根据 实 
地 调查 发 现 , 由 于 近年 来 气候 干旱 及 降雨 减少 ,导致 
流域 西北 地 区 的 典型 草原 已 退化 成 植被 覆盖 度 低 、 
干旱 少 水 的 荒漠 化 草地 或 沙 地 ,使 得 草地 的 平均 ET 
强度 整体 偏 低 。 


表 3 不 同 土地 利用 类 型 ET 估算 结果 
Tab.3 Estimated ET results of different land 


use types /(mm*d^!) 
ds 城乡 工矿 ”未 开发 
平均 值 1.38 1.89 1.11 4.43 1.12 1.32 


4 结论 


针对 遥感 能 量 平 衡 方法 佑 算 区 域 ET. 值 的 空间 
BOSC PE Tey el, LE METRIC 模型 为 例 ,基于 7 与 LAT 
FAIA, FAR T, 频率 直方 图 占 比 前 5% 的 最 高 
温度 区 LAT 频率 直方 图 占 比 后 596 的 裸 土 区 为 干 
OT, 频率 直方 图 占 比 后 5% 的 最 低温 度 区 LAT 


元 ,进而 实现 了 区 域 ET 的 估算 。 经 涡 度 相关 系统 
实测 水 汽 通 量 数据 验证 ,校准 期 和 验证 期 的 区 域 ET 
估算 精度 平均 为 15% ,与 以 往 研究 成 果 相 比 , 结 

相对 合理 ,提取 方案 在 一 定 程度 上 可 为 METRIC 模 
型 主观 确定 干 湿 限 提供 帮助 ;将 ET 按 土地 利用 分 
类 提取 , 研究 区 不 同 土地 利用 类 型 为 : ET a. > 


用 了 林地 > ET yw > 未 开发 和 用 土地 > 上 了 gz 工矿 居民 用 地 ^ 
ET uno AROS DOE ER ET 精准 估算 研究 具有 一 

结合 本 文 研究 结论 ,对 METRIC 模型 估算 区 域 
ET 面临 的 和 干 湿 限 提取 问题 有 如 下 认识 和 思考 :四 
植被 覆盖 良好 、 水 分 充足 的 区 域 或 开放 的 水 体 是 潜 
在 的 “ 湿 点 ”,“ 湿 点 ” 像 元 的 地 表 净 辐射 量 与 土壤 热 
通 量 差 值 几 乎 全 部 转化 成 了 潜 热 通 量 , 使 用 此 类 方 
法 的 前 提 是 ,研究 区 必须 存在 符合 “ 湿 点 ”要 求 或 标 
准 的 像 元 。@ 目前 研究 者 受 认 知 背景 和 下 垫 面 条 
件 的 影响 , 干 湿 点 选取 标准 尚未 统一 , 易 出 现 因 空间 
歧义 导致 的 87 估算 误差 。 本 文 提取 模型 所 需 的 干 
湿 限 是 基于 7. 和 LAI 相关 趋势 图 ,但 估算 的 ET 精 
度 仍 有 提高 的 潜力 。 笔 者 认为 ,深入 研究 干 湿 限 内 
涵 和 标准 ,并 讨论 影响 干 湿 限 分 布 的 主 控 因 子 ,能 对 
模型 的 发 展 起 到 帮助 ,这 也 是 我 们 下 一 步 研究 工作 
的 重点 。 


国家 野外 科学 观测 研究 站 陈 世 蔷 研 究 员 在 涡 度 
统 水 汽 通 量 数据 等 方面 提供 的 帮助 。 
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Dry-Wet Edge Based on Land Surface Temperature and Leaf Area 
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Abstract: The METRIC model was chosen to reveal the spatial ambiguity issue encountered during the calculation 
of regional evapotranspiration (ET) by remote sensing energy balance method. The study was based on the trend 
graph of surface temperature ( T,) and leaf area index (LAI). The dry pixel was chosen as the bare soil area where 
the T, frequency histogram accounted for the top 5% of the high-temperature region, and the LAI frequency histo- 
gram accounted for 5% of the region. The wet pixel was chosen as the wetland or farmland where the T, frequency 
histogram accounted for the last 5% of the low-temperature region, and the LAI frequency histogram accounted for 
the top 5% of the region. After that, the measured data of the eddy covariance system were used to verify the esti- 
mation accuracy of the model. The results showed that the average estimation accuracy of evapotranspiration in the 
calibration period and verification period was 15% , which was relatively reasonable compared with the previous re- 
search results. The selection scheme could help the Mapping Evapotranspiration with Internalized Calibration 
(METRIC) model to subjectively determine the dry-wet edge. The evapotranspiration values of different land use 
> ET > ET. > ET, 


forest arable land unused land 


types were calculated, and they were in an order of ET 


ET, ET 


urban and rural industrial land grassland 9 


> 
which could be referred in researching the precision estimation of regional remote 
sensing evapotranspiration. 
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